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ABSTRACT 

The white-light continuum emission of a solar flare remains a puzzle as regards its height of formation 
and its emission mechanism(s). This continuum, and its extension into the near UV, contain the bulk 
of the energy radiated by a flare, and so its explanation is a high priority. We describe a method 
to determine the optical depth of the emitting layer and apply it to the well-studied flare of 2002 
July 15, making use of MDI pseudo-continuum intensity images. We find the optical depth of the 
visible continuum in all flare images, including an impulsive ribbon structure to be small, consistent 
with the observation of Balmer and Paschen edges in other events. 
Subject headings: Sun: flares — Sun: photosphere 



1. INTRODUCTION 

The first-observed feature of a solar fiare, the "white 
light" continuum, remains enigmatic to the present day. 
Nevertheless it (and the related near-UV continuum) 
contain the bulk of a flare's radiant energy, so that an 
understanding of how it forms would help a great deal in 
our understanding of flare physics. Indeed, understand- 
ing the most important component of the energy of a 
fiare must in fact be the single m ost important problem 
I in fiare physics (e.g., lNeidiglll989l ). 

The difficulties limiting our observational knowledge 
are several. First, the continuum excess over the bright 
photospheric emission only rises to some tens of percent 
even in the most powerful flares, and there are both prac- 
tical (the image contrasts of sunspots) and intrinsic (con- 
vection and other solar "noise," such as the p-modes) lim- 
its to the photometric precision. Also the strong variabil- 
ity both in space and in time of the white-light emission 
has made it difficult to obtain a good characterization 
of the spectrum of the emission. The data that do exist 
suggest two types of white-light flare continuum: those 
associated with the impulsive phase and which correlate 
with particle signatures, and more gradual brightenings 
which have featureless continuum (Boyer et al., 19850. 
We do know that the white-light continuum can 
I have a strong association with the hard X-ray emis- 
, sion in t he impulsive phase of the flare, both 
temporal (iRust fc Hegwerl 119751 : iHudson et all 119921: 
[Neidig & Kane" "1993") and spatial fe.g., iMetcalf et aU 
|2003; Xu et al. 2006). This indeed was one of the first 
motivations for the "thick target model", which as- 
sociates the chromospheric and photospheric effects of 
a fiare with the energy losses of high-energy particles 
(|Naiita fc Or ralllll970l : [Svestkalll970[ )^nd specifically the 
electrons in the 10-100 keV range (jHudson 1972). The 
energetic s of the ele ctrons matches well (Hud sonI 119721 : 
[Fletcher et al.l l200 7a). at least to the extent that we 
understand the energy in the white-light and UV con- 
tinuum. Recent observations of the bolometric lumi- 
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nosity of a solar fiare (fWoods et al.ll2004 : iKretzschmarl 
120081 : iQuesn el et al.l [201 Of ) have proven to be consistent 
with the idea that the white-light continuum and its 
UV extension dominate the flare luminosity, and fur- 
thermore that this luminosity appears in the impulsive 
phase. Thus mechanisms that can explain white-light 
continuum formation must have a close link to the non- 
thermal electrons responsible for hard X-ray emission in 
the impulsive phase. 

The theoretical quandary that has blocked a full under- 
standing is that the 10-100 keV electrons of the impulsive 
phase cannot penetrate as deeply as the non-fiaring vis- 
ible photosphere (the level for which the optical depth 
at 5000A ^5000 = 1, where one might naively expect the 
white-light flare continuum to ori^nat e). This p r oblem 
was exacerbated in the extreme when IXu et al.1 ()2004l ) 
found significant contrast for fiare emission even at 1.56/i, 
nominally the "opacity minimum" region of the spec- 
trum, which according to standard modeling would form 
actually below rsooo = 1- Several possible solutions to 
this quandary have been proposed, viz: 

1. Over- ionization in the chromosphere: at the stop- 
ping depth of electrons of sufficient total en- 
ergy, excess ionization can enhance the cont i nuuin 
adequately acros s the spectrum (jHudsonl [19721 : 
[Aboudarham & H enoux[ [l986): 

2. Radiative back- warming: hydrogen recombination 
radiation frpm the primary stopping height of the 
thick-target electrons heats the photosphere it- 
self sufficiently to pro duce the observed continuum 
([Machado et al.[[T989[ ): 

3. Proton energy losses: protons at a few MeV energy 
penetrate more deeply (Naiita & Orrall 1970); 

4. Wrong model: the thick-target model does not ap- 
ply, and both the hard X-ray emission and the visi- 
ble continuum come from a deeper layer e xcited by 
some other mechanism (e.g. [Unsold [1968] ). 

Of these possibilities items (1) and (2) make clear pre- 
dictions for the optical depth of the visible continuum; 
an overionized layer in the chromosphere will produce an 
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optically thin spectrum of recombination radiation, while 
backwarming would more closely have a blackbody spec- 
tral distribution. What spectroscopy does exists favors 
(1) to a certain extent, because the best impulsive-phase 
optical spectroscopy suggests the presence of the Balmer 
jump and even the Paschen jump (Neidig & Wibor^ 
[1984). The other categories of explanation (3) and (4) 
are more problematic since so little theoretical work has 
been done, but it is clear that the classical thick-target 
model needs revision because the modern data require 
beam intensities greater than seem physically plausible 
([Fletcher fc Hudsonll2QQ8l : iBrown et al.il2QQ9f ). 

The backwar ming mechanism dese rves special men- 
tion, because as lMachado et all ()1989f ) point out, it must 
happen to a certain extent: Balmer radiation from the 
chromosphere emitted downwards can penetrate freely 
to the upper photosphere (the temperature-minimum re- 
gion), heating the medium up to radiative equilibrium. 
Then the usual H~ opacity at a slightly elevated tem- 
perature (say 300 K) could produce the white-light flare. 
This process would be independent of the optical depth of 
the primary emitter in the chromosphere and, if it were 
small, would permit the Balmer and Paschen edges to 
appear in the spectrum (though diluted). The logic be- 
hind this theory rests upon steady-state ID atmospheric 
modeling, so that if 2D spatial structure on the scale of 
the chromosphere/photosphere height difference exists, 
spatial structure will appear in the backwarmed photo- 
sphere; similarly time variations would introduce discrep- 
ancies between the chromospheric UV source and the 
photospheric response. In a backwarmed photosphere 
there would be no spatial structure on scales smaller than 
that of the photosphere/chromosphere height separation. 

In this paper we address the essential property of op- 
tical depth directly, by analysis of SOHO/MDfl images 
of the well-observed flare of 2002 July 15. The essence of 
the technique is simply to correlate the known intensity 
structure of the photosphere at the location of the flare 
brightening with the brightening itself, and to interpret 
this in terms of a simple slab model for the radiative 
transfer (Boyer et al. 1985). To our knowledge no anal- 
ysis of this kind has been carried out before. 

2. THE X3 FLARE OF 2002 JULY 15 

This flare (GOES class X3) famously showed a multi- 
turn helical structure formed in its plasma ejection, as 
observed in TRACE UV images thou ght to represent 
mainly Civ emission ()Liu et al.l l2003f ). The chromo- 
spheric and photospheric properties of the event were 
recorded with Mees Observatory imaging spectroscopy 
in Ha, and with white-light in iagery from th e Imaging 
Vector Magnetograph at Mees (|Li et al.ll2005f ). Unfortu- 
nately RHESSI was at orbit night during the impulsive 
phase, so there are no hard X-ray images available, but 
Owens Valley microwave data were available. Figure [T] 
arrays the flare white light, soft X-ray, and OVSA 1- 
18 GHz microwave emissions as time series for reference; 
Li et al. have estimated energy fluxes for this compari- 
son. To define the impulsive phase Figure [T] also shows 
the GOES time derivative, since RHESSI hard X-ray 
data were not available. 
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Fig. 1. — Time histories for the flare of 2002 July 15, normahzed 
to the individual maxima. The thin black line is GOES 1-8A, and 
the blue line its time derivative. The red points are white-light 
fluxes, total (upper) and the ribbon-like feature only (lower), both 
scaled to the peak of the total. The orange line is the OVSA 18 GHz 
flux density. 

The data we utilize in this study are the MDI "pseudo- 
continuum" intensity data, actually narrow-band sam- 
ples of the continuum near the magnetically sensitive 
photospheric absorption line of Nil at 6768A (Ding et 
al. 2003). The basic data are 1024x1024 time-averaged 
images at 1-minute cadence in a 10.5' square field of view 
at disk center (pixel size 0.615^'; diffraction limit 1.22A/D 
= 1.36^'). The telescope has excellent pointing stability 
and the noise in a given pixel is predominantly solar in 
origin - broad-band variations from convective motions, 
plus the p-modes. The photometric accuracy of the data 
was further improved by applying the seco ndary flatfields 
detailed in the MDI fiatfield repository (jPotts fc Diverl 
I2009f ) Fl. Figure El shows the fiatfield-corrected images. 

Our analysis focuses initially on the crescent-shaped 
fiare region visible clearly in the 20:30 UT image 
(Frame 2 of Figure [2]). We then extend the technique 
to a generalised case. This part of the 15 July fiare is 
particulary easy to analyse as it has the advantage of 
passing over regions with a wide range of photospheric 
intensities, and it proves simple to construct a low- noise 
photospheric background for it. The continuum emis- 
sion has an elongated structure which might be misin- 
terpreted as a loop in projection, but which the multi- 
wavelength observations of Li et al. (2005) show clearly 
to be the eastern ribbon of the fiare (see their Figure 7, 
in particular the frames showing the Ha blue wing at 
20:03:32 and 20:04:01 UT for reference). The Li et al. 
white-light data show a good match to the initial bright 
microwave source (the orange line in Figure [T]), which 
our less- frequent sampling misses. 

Distributions of the intensity along an arc defining the 
midpoint of this ribbon source, and along the same path 
on the prefiare image (one minute prior) are shown in 
Figure El these will be used together with a simple model 
to argue that the fiare optical depth is small. However 
a first indication that this is the case is simply that the 
photospheric structure is clearly superposed on the fiare 
image. 

3. SIMPLE MODEL 



^ The Michelson-Doppler Imager ([Scherrer et al.iri995 li') on board 
the Solar-Heliophysical Observatory spacecraft. 



"^ Ihttp: //soi . Stanford. edu/sssc/MDI_continuum_hr_flatfields/flatfield! 
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Fig. 2. — The time evolution of the flare at one-minute cadence. The upper images show the flatfleld-corrected MDI continuum images, 
and the lower images show the difference between these and a background photospheric image reconstructed by interpolation, color table 
reversed for clarity. 



Our simple model of the radiative transfer in the white- 
light flare consists of a homogeneous slab (optical thick- 
ness r, source function Si) located above an optically- 
thick photosphere at temperature Tq. This analysis fol- 
lows that of Boyer et al. (1985). We assume that the 
photosphere (in a given pixel) has intensity Iq = 5*0 be- 
fore the flare, where is the source function (the Planck 
function if in LTE). Here we explicitly ignore backwarm- 
ing (see Section [1]), assuming that the photosphere does 
not change significantly during the flare. As will be seen, 
this does not contradict the data we discuss for this flare. 
The observed intensity of a given pixel during the flare 
is /i?, which in this model consists of a combination of 
(attenuated) photospheric emission and direct flare emis- 
sion. The observed brightness during the flare is given 
generally by 



If 



(1) 



where Si is the source function for the flaring emission 
layer. 

We illustrate the analysis for r with a simulated white- 
light flare resembling the part of the event we analyze be- 
low. The model photosphere contains weak modulations 
patterned after granules, and a dark 'pore' region. Con- 
sider the intensity profiles of the model flare brightening 
as shown in Figure IH This shows the spatial cross sec- 
tion of three simple examples, each with the same Gaus- 
sian flare emission profile (upper right panel in red), but 
with different optical depths in the slab and thus differ- 
ent emergent intensities (green). Each 'flare' is arranged 
so that the flare brightening has a 20% increase relative 
to the mean photosphere. The lower graphs show the 
total emission before (blue) and during (green) the flare, 
and the upper graphs (red) show the increase in inten- 
sity caused by the flare. In the case with large optical 
depth (r = 1) the photospheric emission is significantly 
attenuated when observed through the flare, and so ap- 
pears with a lower amplitude superimposed on the flare 
emission. If the photospheric emission is subtracted from 
the total flare emission then the result is a combination 
of the fiare emission with an inverted image of the pho- 
tosphere, having a negative correlation with the origi- 
nal photosphere. The right-hand column shows the case 
where the flare has very small optical depth (r = 0.01). 
In this case the total emission during the flare is ap- 



proximately the sum of the unattenuated photospheric 
emission and the flare emission. When the photospheric 
background is subtracted from the flare, the result is the 
true flare intensity profile, with no inverted component 
from the photosphere. It can be seen from this that if it 
is possible to know the intensity of the photosphere un- 
derlying the flare then the optical depth of the flare can 
be determined directly from the degree of spatial corre- 
lation between the flare brightening and the background 
image. 

This is reinforced by the upper graph, where the pho- 
tospheric intensity is subtracted from the flare frame; 
this is also shown as an image in the lower left panel. If 
the flare were significantly optically thick then the varia- 
tions on the attenuated photosphere viewed through the 
flare would be smaller than those viewed directly, and 
therefor this difference image would contain an inverted 
represention of the photosphere, as can be seen in the 
middle column of Figure IH The absence of this tells us 
that the optical depth r must be small, and furthermore 
that the background photosphere has not varied drasti- 
cally during this minute. 

The analysis aims at determining the quantities r and 
Si = bB(Ti) at each pixel of the continuum brighten- 
ing; b is the departure coefficient and B(Ti) the Planck 
function at the source layer temperature Ti. In gen- 
eral, as suming -^"^-^ ( local thermodynamic equilibrium; 
see e.g. lMihalaslll978f ). the opacity and source function 
are related by Kirchoff's law I = /t:S. If LTE held and 
b = 1 we could deduce the two unknowns for each pixel 
directly. In practice we cannot assume LTE, unfortu- 
nately; indeed the phenomenon of white-light continuum 
emission in a flare requires a highly atypical application 
of radiative-transfer theory to the physics of the solar 
atmosphere because of the extreme conditions. The re- 
laxation of the LTE condition for the continuum would 
mean that an additional unknown equivalent to the non- 
LTE departure coefficient must be introduced. We deal 
with this here by using the heuristic approximation de- 
scribed in the following section. 



4. ANALYSIS 



4.1. 



Heuristic opacity model 

We note that an idealized backwarming model for the 
flare brightening would also match this observation, but 
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Fig. 3. — Intensity profile along the selected part of the flare, based on Image Frame 2 (20:04:30 UT). Left upper, the flare brightening 
excess; left lower, the photospheric structure underlying the emission; right, traces along the structure. Here green shows the flare 
brightening, blue the reconstructed photosphere, and red the flare excess. 



Flare continuum opacity 



5 




Fig. 4. — Simulated white-light flare emission proflles compared 
to the underlying photospheric emission for different values of the 
optical depth (r). The variations on the photosphere represent 
granules and a small sunspot, and all intensities are normalised to 
the mean photospheric intensity. Color coding is the same as in 
Figure [3] 
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Fig. 5. — Left: Calculated flare intensities for a = 0.004 and 
a = 0.25 for the path shown in Figure [3] after passing through a 
high-pass fllter with cutoff of 3.5 Mm FWHM. The bold line sec- 
tions are the regions selected for the correlation calculations where 
signiflcant flare brightening is combined with large photospheric 
variation. Notice the correlation between the derived flare inten- 
sity and the photosphere intensity for the a = 0.25 case. Right: 
Correlation of the derived flare intensity with the underlying pho- 
tospheric emission as a function of the opacity parameter a. The 
correlation is zero for a = 0.0037, and the lower panel shows the 
optical depth of the derived flare for this value of a. 

we can rule this ou t because of the fine structure present 
in the fiare image. Ilsobe et"aLi (2007) discuss this aspect 
of white-light flare geometry in detail. In our model, 
the height of the slab is comparable to the width of 
the ribbon-like flare emission, which varies substantially 
along its length and is unresolved in its narrowest parts. 
Furthermore, the flare region we discuss has a lifetime 
less th an 60 s, whereas the cooling time in the VAL_C 
model (|Vernazza et al.lflQSlI ) is about 80 s at the pho- 
tosphere; we calculate the timescale as eH/Cq, where £, 
H, and Cq are the thermal energy density, scale height 
of the photosphere, and photospheric luminosity. This is 
not a decisive argument because of the low time resolu- 
tion of our observations. 

In the case that r is small, we can simplify equation [T] 
by taking a linear expansion of the exponential function 
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If- Soil- r)^ Sir 



(2) 



This approximation gives errors of less than 4% for r < 
0.25, and less than 0.5% for r < 0.1. 
As both r and are unknown Equation [2] cannot 



be solved from just the knowledge of If and So- We 
therefore need to postulate an heuristic relationship be- 
tween r and Si , which would be related to the physics of 
the emission and absorption in the slab region to make 
the problem tractable. If the physical situation in the 
sourc e region (the slab) is characterized by overioniza- 
tion ()Hudsoni ri972) and heating, we would expect r and 
Si to correlate with each other, although the overioniza- 
tion would tend to increase the departure coeflicient for 
the continuum, and hold the temperature Ti to a lower 
level. The simplest relationship would be r = const., but 
this is clearly too restrictive; it would require that faint 
flares perfectly match the photosphere. 

We have therefore adopted r = a{Si — So) as the sim- 
plest functional dependence of opacity on the physical 
conditions in the source. This is consistent with our 
physical argument regarding the effects of non-thermal 
particles heating the flare region. Particle heating would 
generally be expected to increase both the temperature 
and the ionisation faction, leading to strong continuum 
radiation, e.g. in the Balmer free-bound continuum. The 
total energy radiated from the flare volume when the op- 
tical depth is small would then be proportional to the 
physical depth of the heated region and the source func- 
tion in that region. The net radiation which is observed 
will fall between the difference between this and the un- 
derlying photospheric radiation intensity, hence the rela- 
tion above. 

With this heuristic relationship, we flnd a quadratic 
equation for Si and solve it as 



Si = So 



I -So 



(3) 



The value of the constant a, which sets the optical depth, 
is unknown, but we can determine it by considering the 
spatial distribution of the flare emission. Note that a is 
not the optical depth itself, but a parameter describing 
its physics. As can been seen from Figure [3l the photo- 
sphere has more structures at small scale than the flare 
brightening. Recall also our assertion that the flare in- 
tensity is not correlated with the underlying photospheric 
intensity before the flare. The flare source function 
can be calculated using Equation [T] for a range of values 
of a and hence r. If the variation of this derived flare 
source function at small scales is compared to the inten- 
sity of the underlying photosphere interpolated between 
preceding and following image frames, the two should be 
uncorrelated (except for random coincidence) when the 
correct optical depth is used in the calculation. 

4.2. Estimation of the opacity parameter a 

The key to the method described above it to ob- 
serve the attenuation of the photospheric variation when 
viewed through the flare. In order to do this we need to 
identify regions where the following are true: 

1. We must be able to determine the underlying pho- 
tospheric emission at the time of the flare. This 
requires reconstructing the photosphere by inter- 
polating surrounding non-flare image frames. The 
error on this can be estimated by comparing the 
reconstruction to the observed photosphere in re- 
gions outside the flare. 
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Fig. 6. — Selection of regions for analysis. Left: Map of the photospheric RMS variability; middle: difference image showing the flare 
brightening, with contours showing the selected regions; right: the union of these selections for analysis. The label "2" identifles Image 
Frame 2 at 20:04:30 UT. 



2. There must be significant spatial variation in the 
photospheric intensity over the flare region, which 
is much larger than the error on our photospheric 
reconstruction. 

3. There must be significant flare brightening, again, 
much larger than the noise on the photospheric re- 
construction. 

4. The variation of the photospheric intensity needs 
to be large as possible in comparison to the co- 
spatial intensity variation of the flare at the scales 
considered. This is necessary to reduce the effects 
of chance correlations 

The flare cross section shown in Figure [3] is a partic- 
ularly good example that meets these criteria. In this 
case the reconstructed photospheric background was par- 
ticularly accurate as it could be generated by a linear 
interpolation between image frames 3 and 6 (see figure 
[2]) due to the rapid evolution of the flare brightening. 
The error on this reconstruction was ~0.5% of the mean 
photospheric intensity, estimated from the scatter of the 
reconstructed photosphere in regions away from the flare 
to observations. The profiles through this ribbon are 
shown in the right-hand panel of Figure [3l The photo- 
spheric variations are generally at a smaller scale than 
the flare brightening, so to exploit this we passed the 
data though a high-pass filter with cut-off of 3.5 Mm 
FWHM, and chose the areas of high photospheric vari- 
ation manually. These are shown by the bold portions 
of the filtered cross section shown in Figure [51 A flare 
model profile was then generated using the reconstructed 
So and Equation [3] for a range of different values of the 
opacity parameter a, and the correlation between this 
and the photospheric background was calculated. Fig- 
ure [5]shows the results, with reconstructions for a = 0.25 
and a = 0.004 shown in the left panels, and the variation 
of the correlation with a shown on the right hand panel. 
For the case a = 0.25 the generated flare profile shows a 
significant positive correlation with the underlying pho- 
tosphere, which completely disappears when when a is 
small. According to the arguments in Section [3l this im- 
plies small optical depths. In this case the correlation 
goes to zero when a = 0.0036. From our assumption 
that r = a{Si — So) this gives an average optical depth 



TABLE 1 

Summary of the results from different image frames 



Frame 


ao 


Mean r 


Photospheric 
rms error 


AS/S (%)« 


2 


0.0037 


0.015 


0.78 


8.5 


3 


0.0131 


0.033 


1.17 


11.5 


4 


0.0309 


0.042 


1.28 


7.2 


5 


0.0119 


0.024 


1.31 


5.9 


6 


0.0233 


0.032 


1.35 


5.2 


7 


0.0066 


0.016 


1.32 


4.7 


all frames 




0.028 







"AS/S is the flare excess brightness, from the selected regions, 
normalized to the non-flaring photosphere. 

for the flare of r = 0.0016, with the variation of r along 
the path shown in the the lower panel. This results from 
the low photospheric noise. 

The above example indicates very small optical depth, 
but since the method assumes no correlation between 
the detailed photospheric intensity and flare brighten- 
ing, random correlations could also give that result. To 
test this we need to look at a larger data set. To auto- 
mate this algorithms were constructed that isolated re- 
gions where the four criteria listed above were true. An 
example of this from Frame 2 is shown in Figure [6l This 
process was carried out for all image frames and the re- 
sults are summarised in Table [TJ As can be seen, the 
mean optical depth of the flares in each frame is small, 
with a maximum of 0.042. Notice that the value of r is 
not related to the average flare brightening, and is low- 
est for the frame which has the lowest photospheric error 
(Image Frame 2). 

4.3. Error analysis 

The major error source in the above procedure is the 
error in the estimate of the photospheric brightness in 
regions under flare areas. Although the average of this 
error will be zero by assumption, the spurious structure 
it introduces to the photospheric reconstruction gives a 
systematic error to the derived value of a, and hence to 
the derived optical depth, always tending to increase the 
derived value. To see why this should be true consider 
the effects of adding a pattern of random noise to the re- 
construction of the non-flare photosphere. The observed 
intensity of the solar surface with the flare present will 
not contain this error term, so it will appear that the 
flare is masking it, suggesting that it is more optically 
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Fig. 7. — The effect of noise in the photospheric reconstruction 
on the derived value for a. Extrapolating this graph back to the y 
axis allows the zero noise case to be estimated. 

thick. 

In order to evaluate the significance of the photospheric 
noise, we add Gaussian random noise of various (known) 
amplitude levels to the reconstructed photosphere. First 
the RMS error on the photospheric reconstructions was 
measured for each image frame by comparing the recon- 
structed photospheric image in a region away from the 
flare emission to the observed photosphere; results are 
shown in Table [TJ Normally distributed noise with a 
range of known standard deviations was then added to 
the reconstructed photosphere and the calculation of a 
repeated. This was done multiple times with different 
random noise for each frame and each value of the noise 
amplitude, in order to avoid random correlations between 
added noise and the photosphere. The variation in the 
derived value of a versus the total photospheric error for 
each frame is shown in Figure [71 It can be seen that 
as the error in the photospheric noise increases, the de- 
rived value of a also increases as expected. Note that the 
values of a determined in image frames where the flare 
brightness is larger (see Table [T]) are much less sensitive 
to the effect of the noise. 

Using this graph it is possible to extrapolate back to a 
'zero error' value for a. If this is done for each frame us- 
ing a simple polynomial extrapolation we get an average 
value for all the data of a = —0.0001 ± 0.01, so to within 
the accuracy of these measurements the optical depth is 
effectively zero. 

4.4. Discussion 

Our direct analysis of the flare emission from seven im- 
age frames, comprising a total of 1200 pixels where the 
flare emission significantly enhanced the surface bright- 
ness, has given us the result that the optical depth is 
0.028 ±0.01. This value however should be regarded as 
an upper limit on the opacity of the flaring regions; when 
the effect of photospheric noise is considered the optical 
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depth becomes too small to measure, and certainly less 
than 0.01. 

As a result of this the assumption of our heuristic opac- 
ity model that r is depends on the flare source function 
becomes unimportant, as the optical depth is so close to 
zero that a more appropriate model for the emission be- 
comes /i? = /So + /i, where Ii is the flare emission. The 
flare excess simply adds to the photospheric emission. 

5. CONCLUSIONS 

The optical depth of the white-light flare regions we 
have studied is very small, indistinguishable from zero 
in this study, and in any case less than ~0.01. We infer 
from this that the flare must be of low density and hot, 
almost certainly far from LTE. The temperature can- 
not be determined from this but is generally constrained 
by the new flare bolometric observations ([W oods et al.l 
120041 : lKretzsch mar|[2QQ8l : IQuesnel eFaIll201Qft : many au- 
thors suggest a value near 10^ K (e.g., Hudson et al. 
2010). These results definitely tend to reduce the im- 
portance of photospheric backwarming in our under- 
standing of white-light flare emission. First, the low 
density and high temperature imply an emission source 
high in the atmosphere, consistent with stopping depths 
of low-energy electrons but inconsistent with stopping 
depths of the electron s required for backwarming models 
([Fletcher et al.|[2"007bf ). Second, the flare image scales (in 
this case limited by MDI resolution) are of order 1 Mm. 
If backwarming contributed significantly to the emission 
then the maximum height of the emission must be around 
half of the feature size (simple geometric ray model), 
which would be too low for an optically thin case (see 
llsobe et al.ll200 7). 

Because the MDI data represent averages over one- 
minute intervals, and because they represent only a nar- 
row sHce of the true continuum, this result should be 
considered as a preliminary one. Newer data with better 
image cadence and spatial resolution (Hinode or SDO in 
space, or a variety of ground-based instruments leading 
up to ATST) should be applied to this interesting prob- 
lem. If confirmed, this result suggests that the generally 
accepted picture of flare energy storage in the corona, 
with flare effects in the lower atmosphere derived from 
this energy reservoir, must be correct. Note that this is 
the usual assumption, but that it has not been easy to 
establish observationally. 
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